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ABSTRACT 

We study the spectral evolution of 13 short duration Gamma Ray Bursts (GRBs) detected by 
the Gamma Burst Monitor (GBM) on board Fermi. We study spectra resolved in time at the 
level of 2 - 512 ms in the 8 keV-35 MeV energy range. We find a strong correlation between 
the observed peak energy -Epeak and the flux P within individual short GRBs. The slope of 
the -Epcak oc P s correlation for individual bursts ranges between ^0.4 and ~1. There is no 
correlation between the low energy spectral index and the peak energy or the flux. Our results 
show that in our 13 short GRBs -E pGa k evolves in time tracking the flux. This behavior is 
similar to what found in the population of long GRBs and it is in agreement with the evidence 
that long GRBs and (the still few) short GRBs with measured redshifts follow the same rest 
frame -E pca k — ^iso correlation. Its origin is most likely to be found in the radiative mechanism 
that has to be the same in both classes of GRBs. 

Key words: Gamma-ray: bursts — Radiation mechanisms: non thermal 
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1 INTRODUCTION 

Short Gamma Ray Bursts have been a challenge since the finding of 
their spectral diversity with respect to the class of long GRBs (e.g. 
see Nakar 2007; Lee & Ramirez-Ruiz 2007 for recent reviews). 
Short GRBs have optical and X-ray afterglows, and in few cases 
they show also X-ray flares, similar to those of long bursts but 
scaled by their fluence (e.g. Gehrels et al. 2008; Nysewander et 
al. 2008). However, the class of short GRBs is somewhat heteroge- 
neous for what concerns the prompt emission properties (e.g. there 
are short GRBs followed by a faint extended emission - Norris & 
Bonnell 2006; Donaghy et al. 2006; Norris & Gehrels 2008) or the 
host galaxies properties (short GRBs are found in almost all galaxy 
types - e.g. Berger 2009). Recently, this picture has been also com- 
plicated by the detection of long GRBs at very high redshifts like 
GRB 080913 atz = 6.7 (Greiner et al. 2009), and GRB 090423 at 

2 ~ 8.2 (Salvaterra et al. 2009, Tanvir et al. 2009) which have an 
intrinsic duration of less than 2 s. Classification schemes of short 
GRBs that try to merge all these evidences have been proposed 
(Zhang et al. 2007). 

For what concerns the prompt emission, it was discovered 
through the BATSE sample that short GRBs are spectrally harder 
than long bursts (Kouveliotou et al. 1993). Detailed analysis of 
the time integrated spectra of short GRBs (Paciesas et al. 2003, 
Ghirlanda, Ghisellini & Celotti 2004, Ghirlanda et al. 2009) have 
also shown that this spectral difference is due to a harder low energy 
spectral component in short GRBs (but see Nava et al. 2010). These 
studies also revealed that, on average, the first 2 seconds of emis- 
sion of long GRBs have similar spectral properties of short bursts 
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(as also found from the comparison of the variability patterns - 
Nakar & Piran 2002), thus suggesting the presence of a common 
emission mechanism that operates in these two kind of sources. 

Considered individually, long GRBs show a strong spectral 
evolution with two possible behaviors: the peak energy Speak of 
the vF v spectrum decays in time ("hard to soft" evolution) or fol- 
lows the variation of the flux ("tracking" evolution; e.g. Band et 
al. 1993; Ford et al. 1995). This emerged from the time resolved 
spectral analysis of long GRBs observed by BATSE. Recently Lu, 
Hou & Liang (2010) reported that the hard to soft spectral evolu- 
tion is dominant in long bursts (at least initially, for 2/3 of their 
sample of 22 single pulses BATSE GRBs). To date, however, no 
detailed study of the possible spectral evolution of short GRBs was 
performed. 

Another recent issue concerning the prompt emission of GRBs 
is the nature of the spectral-energy correlations discovered for the 
sample of long GRBs with measured redshifts. Amati et al. (2002) 
found that the rest frame peak energy correlates with the isotropic 
energy (the -B pca k — -Eiso correlation) while Yonetoku et al. (2004) 
showed that a similar correlation exists between -E pca k and the 
isotropic luminosity (the i? pca k — ^iso correlation). The open de- 
bate is if these correlations have a physical foundation (Ghirlanda 
et al. 2005; Bosnjak et al. 2008; Ghirlanda et al. 2008; Nava et 
al. 2008; Krimm et al. 2009; Amati et al. 2009) or if they are the 
result of selection effects (Nakar & Piran 2005; Band & Preece 
2005; Butler et al. 2007; Butler et al. 2009; Shahmoradi & Ne- 
miroff 2009). Both the i?p C ak — Ei so and -Bpcak — £iso correlations 
are defined by considering the time integrated spectral properties 
of long GRBs with measured redshifts. Recently, Ghirlanda et al. 
(2009) showed that also short GRBs with measured redshifts fol- 
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low the same -E pea k — Liso correlation defined by long bursts, but 
not the corresponding T pca k — Ei so correlation. 

One of the most stringent result supporting a physical origin 
of the -Epcak — ^iso correlation is that, in long GRBs with mea- 
sured redshifts detected by Fermi, there is a strong Tp Ca k — Li so 
correlation within individual bursts (Ghirlanda et al. 2010) which 
cannot be due to selection effects but must have a physical ori- 
gin. A similar result was reached by Firmani et al. 2009 based on 
the time resolves spectral analysis of long GRBs detected by Swift. 
Several interpretations for the Speak — Li so (or -E pca k — -Eiso) corre- 
lations have been proposed (e.g. Yamazaki et al. 2004; Lamb et al. 
2005; Rees & Meszaros 2005; Levinson & Eichler 2005; Toma et 
al. 2005; Eichler & Levinson 2005; Barbiellini et al. 2006; Thomp- 
son 2006; Ryde et al. 2006; Giannios & Spruit 2007; Thompson 
et al. 2007; Guida et al. 2008; Panaitescu 2009), but there is no 
general consensus about a prevalent idea. 

Therefore, the two topics we aim to address in the present 
work are (1) whether, and how, the spectra of short GRBs evolve 
in time, (2) if also in short GRBs there is a spectral-energy corre- 
lation between the peak energy and flux similar to what found in 
long events. 

Time resolved spectral studies of short GRBs are hampered by 
their lower fluence and duration with respect to long GRBs. More- 
over, time resolved spectroscopy requires an instrument with the 
largest energy coverage and good spectral resolution in the keV- 
MeV energy range in order to constrain the spectral parameters 
of individual (time-resolved) spectra and to follow their tempo- 
ral variation. Swift detected a large number of short GRBs, but it 
is not suited to this aim because of the narrow energy range of the 
BAT instrument (15-150 keV). The Gamma Burst Monitor (GBM; 
8 keW-AO MeV) on-board the Fermi satellite, instead, represents a 
unique opportunity to study, in details for the first time, how the 
spectrum in short GRBs evolves with time. 

In §2 we present the sample of short Fermi GRBs selected 
for the time resolved spectral analysis. The details of the spectral 
analysis are described in §3 and the results in §4. We draw our 
conclusions in §5. 



2 THE SAMPLE 

We consider the 237 GRBs detected by the Fermi GBM up to May 
2010 and whose detection has been published through GCN com- 
munications. In this sample there are 37 short GRBs with observed 
duration (as reported in the GCNs) 5C2s. For the time resolved spec- 
tral analysis we consider the short GRBs with the largest fluence 
and peak flux. This ensures to have, for each GRB, a set of time 
resolved spectra with enough signal to constrain the spectral pa- 
rameters. We select the 14 short GRBs with a fluence larger than 
8 x 10 ~ 7 erg cm -2 (integrated in the 8 keV - 1 MeV energy range) 
and a peak photon flux P pC ak H ph cm" 2 s" 1 . 

We anticipate that for one burst (GRB 100223) the detector 
response files are not present in the archive yet so that its spectral 
analysis was not possible. GRB 081 1 13 has a time integrated spec- 
trum which is consistent with a single power law (von Kienlin et al. 
2008a) and similarly its time resolved spectra are all consistent with 
a single power law. All the other GRBs of our sample have time in- 
tegrated spectra (as reported in the GCN) with a peak in vF v . We 
have included GRB 090308 and GRB 081107 in the sample al- 
though their duration (Tgo) is somewhat longer than 2 s (i.e. 2.2 s 
and 2.11 s respectively), since they have a fluence and peak flux 
larger than our threshold and considering that the division of short 



Table 1. Short GRB sample. The last three digits in parenthesis are the frac- 
tional trigger number assigned to the GRB by the GBM archive. The fluence 
.F_6 is in units of 10 — 6 erg cm -2 . References for the duration (Tgo), flu- 
ence F-a and peak flux P pca k : (1) Guiriec et al. 2009; (2) Guiriec et al. 
2009a; (3) von Kienlin et al. 2009; (4) McBreen et al., 2008; (5) Wilson- 
Hodge et al., 2009; (6) Bissaldi et al., 2009; (7) Goldstein et al 2009; (8) 
Bissaldi et al., 2008; (9) von Kienlin et al., 2008; (10) Bissaldi et al., 2009a; 
(11) Goldstein et al 2009a; (12) von Kienlin et al., 2010; (13) Goldstein 
et al., 2009b; (14) von Kienlin et al., 2008a. Col. 6 lists the range of time 
resolution in milliseconds used to perform the spectral analysis. In the case 
of GRB 090228 the time resolution is 2 ms for the first peak and 32 ms 
for the (fainter) second one of its lightcurve. * GRBs with a precursor. 
GRB 090308 and GRB 081107 are shown separately since they have a 
slightly larger duration than the canonical 2s dividing line between short 
and long GRBs but still have a fluence and peak flux above our selection 
thresholds. 
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and long GRBs at 2 s is not sharp, since both populations have a 
gaussian Tgo distribution extending below and above this time cut. 
These two bursts are reported separately in Tab. Q] Therefore, the 
sample of analyzed short GRBs is composed by 14 events. 

In Tab.Q]we report the names, duration, fluence and peak flux 
of the selected short GRBs (with the corresponding reference - Col. 
4). The redshift is not measured for all the bursts in our sample 
except for GRB 090510 at z = 0.903 (Rau et al. 2009). 



3 SPECTRAL ANALYSIS 

The GBM (Meegan et al. 2009) comprises 12 thallium sodium io- 
dide [Nal(Tl)] and two bismuth germanate (BGO) scintillation de- 
tectors which cover the energy ranges ~8 keV-1 MeV and ~300 
keV-40 MeV, respectively. The GBM acquires three types of data 
suited for spectroscopic studies (Meegan et al. 2009): the "CTIME" 
data consisting of a sequence spectra (binned in 8 energy chan- 
nels) with a time resolution between 0.256 and 1.024 seconds, the 
"CSPEC" data containing a sequence of 128 energy channel spectra 
binned in time with a variable resolution between 1.024 s and 4.096 
s and the "TTE" event data files containing individual photons with 
time and energy tags. 

To the aim of studying short GRBs, TTE data are ideal because 
they allow to choose the temporal resolution according to the burst 
duration and flux. Due to the limited buffer size, TTE data only 
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Figure 1. Correlation between the peak energy E pcai[ and the flux of the 163 time resolved spectra of the 13 short GRBs analyzed in this work (Tab.[T}. 
Different symbols/colors correspond to different bursts (as shown in the legend). The solid line (dashed lines) is the "Yonetoku" relation (and its 3<r scatter) 
of long GRBs with measured redshifts (adapted from Ghirlanda et al. 2009) transformed in the observer frame assuming 2 = 1. The grey (solid and dashed 
lines) are for z = 0.5. 



covers the interval between ~30 s before and ~300 s after the burst 
trigger time. This time interval fully encompasses the duration of 
short GRBs and allows to fit the background spectrum by selecting 
time intervals before and after the burst. 

For the time resolved spectral analysis we used the recently 
released software RMFI10(v33pr7). In order to model the back- 
ground spectrum for the time resolved spectroscopic analysis, we 
selected two time intervals before and after the burst. The sequence 
of background spectra in the two selected intervals were fitted with 
a first order polynomial to account for the possible time variation of 
the background spectrum. Then the background spectrum was ex- 
trapolated to the time intervals selected for the time resolved spec- 
troscopy of individual GRBs. 

For each burst, we jointly fitted the spectra from the Nal de- 
tectors which had the largest illumination by the GRB. The Nal 
detectors were selected for having an angle of position of the GRB 
with respect to the detector normal lower than ~80 degrees. Among 
the two BGO detectors we selected that with the smallest angle 
to the GRB. Four Nal detectors and one BGO were analyzed for 
GRB 100206, 090902, 081213, 090108, 090305, 090308, 081216, 
090228, 090227B; three Nal and one BGO for GRB 090206, 
090328 and three Nal for GRB 081 107. For GRB 090510 we used 
5 Nal and both BGO data. The inclusion of the BGO data extends 
the spectral coverage of the Nal detectors from 1 MeV to ~35 MeV. 



1 http://fermi.gsfc.nasa.gov/ssc/data/analysis/user/ 



This energy extension represents an unprecedented opportunity for 
the spectral analysis of short Fermi GRBs whose time integrated 
spectrum has typically a peak energy of the order of 1 MeV (e.g. 
Navaetal. 2010). 

For each burst we performed a time resolved spectral analysis 
by opportunely changing the time resolution (starting with 64 ms 
time resolution of the TTE data) and we checked that each time 
resolved spectrum had well constrained spectral parameters (i.e. 
the normalization, E pca ^ and a are required to have relative errors 
less than 100%). For the brightest GRBs of our sample (i.e. GRB 
090227B, 090228) we performed a time resolved spectral analysis 
down to the 8 ms and 2 ms timescale. In some cases, the end of the 
slow decline of the light curves did not ensure enough signal to ob- 
tain a good fit so that a coarser time resolution was applied (e.g. a 
single bin of 0.5 s at the end of GRB 090902 was analyzed). In Tab. 
Q] (last column) we report the minimum and maximum time reso- 
lution in milliseconds at which the time resolved spectral analysis 
was performed. 

The model adopted is a power law with an exponential cutoff 
whose free parameters are the low energy spectral index a, the peak 
energy i? pea k (i-e the peak of the vF v spectrum), and the normal- 
ization. We also allowed for a variable normalization factor in order 
to fit together the data of the Nal and BGO detectors. For the time 
resolved spectral analysis we fixed the normalization factors to the 
values obtained from the fit of the time integrated spectra of each 
burst. The values of these normalization factors are all consistent 
with 1, with a deviation of 20% at most. The choice of this model 



© 0000 RAS, MNRAS 000, 000-000 



4 G. Ghirlanda et al. 



- SK 09051(1 (prec) 
. 090510 




10" 7 10" 6 10" 5 10" 4 10" 3 

Flux [erg/cm7s] 

Figure 2. Correlation between £? pC ak and the flux P for the 4 GRB 
with the highest fluence. The upper panel shows GRB 090510 and GRB 
090227B. The precursor of GRB 090510 is shown with the green arrow 
symbols since its spectrum is consistent with a single powerlaw. The lower 
panel shows GRB 081216 (with a precursor, green symbols) and GRB 
090228. Color code (data points, solid and dashed lines) as in Fig.[T] 
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Figure 3. £? poa k versus the photon spectral index of the fitted model 
(cutoff-power law). Same symbols and colors as in Fig. [T] The vertical 
lines represent the limits predicted by the standard synchrotron model (dot- 
ted line) and in the case of cooling (dashed line). 



is motivated by the fact that at high energies the response of the 
BGO rapidly decreases for increasing energy so that it is hard, in 
single time resolved spectra, to constrain the possible presence of a 
power law component of e.g. the Band function (Band et al. 1993) 
which is instead typically fitted to the time integrated spectra of 
long GRBs. We tried to fit the time resolved spectra with the Band 
function but in most cases the high energy spectral index ft was 
unconstrained. This is confirmed also by the recent analysis of the 
brightest three short GRBs of our sample (Guiriec et al. 2010) that 
shows that the fit with a Band function results in unconstrained /3 
in most of their time resolved spectra. Therefore, for homogeneity 
and with the aim of comparing the spectral evolution trends in short 
GRBs, we adopted the same spectral model, i.e. a cutoff-power law 
for all the time resolved spectra. 



4 RESULTS 

We searched for possible correlations between the best fit spec- 
tral parameters of the time resolved spectra of our sample of short 
GRBs. We find that there is a strong correlation between the peak 
energy S pea k and the flux P within individual short GRBs. The 
flux P is calculated by integrating the best fit model of each time 
resolved spectrum in the 8 keV-35 MeV energy range. This cor- 
relation is shown in Fig.Q]for the whole sample of 13 GRBs. We 
could extract for them a total of 163 time resolved spectra. Consid- 
ering the 13 short GRBs individually, the S pca k-P s correlation of 
individual GRBs has a slope that ranges between 0.4 and 1.0. The 
Speak — P correlation is in the observer frame (only GRB 090510 
has a known redshift). For comparison with the population of long 
GRBs, we show in Fig.[T]the Speak — Li so ("Yonetoku") correlation 
of long GRBs with measured redshifts (adapted from Ghirlanda et 
al. 2010) assuming a redshift z = 1 and z — 0.5 (solid black and 
grey lines respectively). Also shown are the boundaries represent- 
ing the 3ct scatter of this correlation (dashed lines in Fig. []}. The 



short GRBs spectral evolution trends are consistent with this cor- 
relation (transformed in the observed frame) as already shown for 
the spectral evolution tracks of long GRBs (Ghirlanda et al. 2010). 
However, we note from Fig. Q] that the Speak — P correlation de- 
fined by short GRBs is similar in slope to the S poa k — £iso corre- 
lation transformed in the observer frame (solid grey and black lines 
in Fig.Q} but it has a higher normalization. 

The most sampled events (i.e. those studied on the shortest in- 
tegration timescales - down to 4 ms in the case of the peak of GRB 
090228) are the 4 GRBs with the largest fluence/peak flux (Tab. 
Q}, shown separately in Fig. [2] Noteworthy, two of them also show 
precursor activity which is shown with different symbols (colors) in 
Fig|2] GRBs are known to show, in ~ 15% of cases (Lazzati 2005, 
Burlon et al. 2008), emission preceding the main episode (i.e. pre- 
cursors). Burlon et al. (2009) analyzed a large sample of bright long 
BATSE GRBs with precursor emission and showed that the same 
spectral evolution is present in the precursors and in the main GRB 
events. They also did not find any relation between the spectral in- 
dex and Speak similarly to what shown in Fig.[3] Here we find that, 
although still based only on two cases, also in short GRBs the pre- 
cursors have spectra which are consistent with the general trend of 
the main event. This result, shown for short GRBs in this paper for 
the first time, points towards a common origin of the precursor and 
the main emission. 

Fig. [3] shows that there is no correlation between the low en- 
ergy spectral index a of the cutoff-power law model and the peak 
energy S poa k. Similarly we do not find any correlation between a 
and the flux P. We instead find, as already shown by the analysis of 
time integrated spectra of short GRBs (e.g. Ghirlanda et al. 2009), 
that they are harder than long ones and this makes a large fraction of 
them inconsistent with the "lines of death" of synchrotron emission 
(with no cooling - vertical dotted line in Fig. [3} and makes all of 
them inconsistent with synchrotron emission by cooling electrons 
(dashed vertical line in Fig.[3j- 



© 0000 RAS, MNRAS 000, 000-000 



Fermi Short GRBs 5 



5 CONCLUSIONS 

The time resolved spectra of individual short GRBs evolve in time, 
and their -E pC ak is strongly correlated with their flux. Furthermore, 
the found correlation is very similar to what already found in indi- 
vidual long GRBs. The GBM data for the brightest bursts allowed 
to analyze their spectra even at the 2-8 ms time resolution for part 
of the duration of these GRBs. With respect to typical peak fluxes 
measured on ~1 s time bins, the peak fluxes reached by our short 
GRBs, measured on a fraction of second, are extreme, i.e. of the 
order of severalx 10 -4 erg cm -2 s _1 (see Fig.[TJ. 

We do not know if these flux levels are reached also in long 
GRBs, since a time resolved analysis with this degree of accuracy 
has not yet been done (but time resolved spectra with a coarser 
time resolution for the brightest long bursts had much smaller peak 
fluxes, e.g. Kaneko et al. 2006). The fact that the Speak — P corre- 
lation is similar suggests that long and short GRBs share the same 
emission mechanism for their prompt emission. Furthermore, the 
emission mechanism should not depend on the progenitor, (e.g. 
fireball-funnel interactions) if long and short GRBs do have dif- 
ferent progenitors. 

For the brightest bursts, the dynamic range of the Speak — P 
is very large, being more than 2 orders of magnitude in flux and 
more than one in S pC ak. Even precursors obey the same S pC ak — P 
correlation, although only two short GRBs in our sample show a 
precursor. 

The three brightest short bursts in our sample reach values of 
Speak significantly larger than those of long GRBs (e.g. Ghirlanda 
et al. 2010). On the other hand, they also reach significantly larger 
fluxes while remaining on the same -Z5 pca k — P correlation defined 
by long bursts. For the same fluxes they do have the same Speak- 

The iJp C ak — P correlation, as well as the analogous Yonetoku 
relation for time integrated spectra of different GRBs, has not yet 
received a convincing and broadly accepted explanation, as men- 
tioned in the introduction. The fact that also short bursts obey it 
makes the search for a convincing interpretation even more com- 
pelling. 
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